Measurement of the Positive Muon Anomalous Magnetic Moment to 0.7 ppm 
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A higher precision measurement of the anomalous g value, a M = (g — 2)/2, for the positive muon 
has been made at the Brookhaven Alternating Gradient Synchrotron, based on data collected in the 
year 2000. The result = 11 659 204(7) (5) x 10 -10 (0.7 ppm) is in good agreement with previous 
measurements and has an error about one half that of the combined previous data. The present 
world average experimental value is a M (exp) = 11659 203(8) x 10~ 10 (0.7ppm). 

PACS numbers: 13.40.Em, 12.15.Lk, 14.60.Ef 



The study of magnetic moments has played an im- 
portant role in our understanding of sub-atomic physics. 
Precision measurements of the electron anomalous mag- 
netic moment, together with those of the hyperfme struc- 
ture of hydrogen and the Lamb shift, stimulated the de- 
velopment of modern quantum electrodynamics and have 
since provided stringent tests of this theory. In this Let- 
ter we report a new measurement of the anomalous mag- 
netic moment of the positive muon, a M = (g — 2)/2, with 
a relative precision of 0.7 parts per million (ppm), nearly 
two times better than our previous work jl, |[ ||. This 
measurement comes from data collected in the year 2000. 
At this level, a M is sensitive to QED, weak, and hadronic 
virtual loops and provides an important constraint on 
extensions to the standard model. 

The principle of the experiment and previous results 
have been given in earlier publications ||. Also, 
detailed descriptions of the (g — 2) superconducting in- 
flector magnet, storage ring magnet, fast kicker, NMR 
system, and calorimeters have been published B . 



The quantity is determined from 



B) 



(1) 



The magnetic field (B) weighted over the muon beam 
distribution is measured by proton NMR. The difference 



frequency uj a between the muon spin precession and or- 
bital angular frequencies is determined by counting the 
number N(t) of decay positrons with energies larger than 
an energy threshold, 

N(t) = N (E)e^ ( ^ [1 + A{E) sm(u a t + ME))} ■ 

(2) 

The normalization Nq, asymmetry A, and phase <f> a vary 
with the chosen threshold. The time dilated lifetime is 
7T Ri 64.4 [is. For muons with 7 = 29.3, the angular dif- 
ference frequency u> a is not affected by electrostatic fo- 
cusing fields in the ring. 

New aspects of the 2000 data taking period include: 
the operation of the AGS with 12 beam bunches, con- 
tributing to a 4-fold increase in data collected as com- 
pared to 1999; a new superconducting inflector magnet, 
which improved the field homogeneity in the muon stor- 
age region; the installation and operation of a sweeper 
magnet in the beamline, which reduced AGS background; 
and additional muon loss detectors, which enable an im- 
proved determination of the time dependence of muon 
losses. Most other experimental aspects of the data tak- 
ing in 2000 were the same as in 1998 and 1999. 

The magnetic field value was obtained from NMR mea- 
surements of the free proton resonance frequency. A trol- 
ley with 17 NMR probes was used to measure the field 
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Multipoles [ppm] 
normal skew 
Quad 0.24 0.29 
Sext -0.53 -1.06 
Octu -0.10 -0.15 
Decu 0.82 0.54 



radial distance [cm] 



FIG. 1: A 2-dimensional multipole expansion of the field av- 
eraged over azimuth from one out of 22 trolley measurements. 
Half ppm contours with respect to a central azimuthal average 
field Bo = 1.451 274 T are shown. The multipole amplitudes 
relative to Bo, are given at the beam aperture, which has a 
radius of 4.5 cm and is indicated by the circle. 



TABLE I: Systematic uncertainties for the ui p analysis 



Source of errors 


Size [ppm] 


Absolute calibration of standard probe 


0.05 


Calibration of trolley probe 


0.15 


Trolley measurements of Bo 


0.10 


Interpolation with fixed probes 


0.10 


Uncertainty from muon distribution 


0.03 


Others 1 


0.10 


Total systematic error on u) p 


0.24 



higher multipoles, trolley temperature and voltage re- 



sponse, eddy currents from the kickers, and time-varying 
stray fields. 

throughout the muon storage ring, typically every three 
days. The trolley probes were calibrated to an accu- 
racy of 0.15 ppm with respect to a standard spherical 
H2O probe, which has an absolute calibration known to 
0.05 ppm. The probes were intercalibrated during the 
data taking period using a single, movable probe plunged 
into the storage region. Interpolation of the magnetic 
field between trolley measurements was based on the 
continual readings of about 150 fixed NMR probes dis- 
tributed around the ring in the top and bottom walls of 
the vacuum chamber. Figure [l] shows a magnetic field 
profile averaged over azimuth. The improved field ho- 
mogeneity obtained with the new inflector relaxes the 
demands on the knowledge of the muon beam distribu- 
tion. 

The field (B) averaged over the muon beam distribu- 
tion in space and time was obtained from two complete 
and largely independent analyses, whose results were 
found to agree to within 0.05 ppm. Its final value is ex- 
pressed in terms of the free proton resonance frequency 
and is given by Wp/(2vr) = 61 791 595(15) Hz (0.2 ppm). 
Table | lists the uncertainties. 

The frequency oj a was obtained from the time distri- 



bution of decay positron counts. The positrons were de- 
tected with 24 lead/scintillating fiber calorimeters on the 
inside of the storage ring, whose photomultiplier signals 
had a typical FWHM of 5 ns. The signals were recorded 
with 400 MHz, 8-bit waveform digitizers (WFD). Pulses 
above a preset energy threshold of about 1 GeV triggered 
the WFD to record at least 16 ADC samples (40 ns). 
Pulses with energies below the hardware threshold were 
recorded if they appeared within the sampling window 
around a trigger pulse. The positron arrival times and 
energies were reconstructed offline from the WFD record- 
ings using two independent implementations of our pulse 
reconstruction algorithm. 

The large positron sample and 30% higher instanta- 
neous count rates than in 1999 required careful consider- 
ation of distortions to the spectrum in Eq. |^. As in the 
analysis of our 1999 data, we have considered (1) positron 
pulses overlapping in time (pileup), (2) coherent betatron 
oscillations, (3) beam debunching, (4) muon losses, and 
(5) detector gain and time instability. 

(1) Pileup of positron signals distorts the time spec- 
trum because of misidentification of the number, ener- 
gies, and times of the positrons. The pileup spectrum can 
be constructed from the data using the extended pulse 
sampling by the WFD described above, and can effec- 
tively be subtracted from the data prior to fitting || . 

(2) As in earlier data taking periods, the weak focusing 
storage ring was operated with a field index n = 0.137, 
well away from beam and spin resonances. The phase 
space for betatron oscillations defined by the acceptance 
of the storage ring was not filled, which resulted in co- 
herent betatron oscillations (CBO) — betatron oscilla- 
tions of the beam as a whole. Since the acceptance of 
a calorimeter varies with the radial muon decay posi- 
tion in the storage ring and the momentum of the pro- 
duced decay positron, the observed positron time and 
energy spectra are modulated with the CBO frequencies. 
The most important modulation is the horizontal one 
with frequency w c bo,h = (l — \A ~ n ) w c = 2ir ■ 466kHz, 
where lo c is the cyclotron frequency. At injection time, 
this modulation affects the terms N , A, and 4> a in Eq. 
at the level of 1%, 0.1%, and lmrad, respectively. 

For n — 0.137, the frequency w c bo,h is approximately 
twice as large as oj a and, hence, the interference frequency 
<^cbo,h — w a is close to u a . Modulations of the asymmetry 
and phase with frequencies w c bo,h - 2 x w„ may, unlike 
the modulation of observed counts, manifest themselves 
as sizable, artificial shifts in the fitted frequency value 
u) a , when not taken into account in the function fitted 
to the data. These shifts are as large as 4 ppm for indi- 
vidual calorimeter time spectra, and mostly cancel in the 
summed spectra owing to the circular symmetry of the 
experiment design. 

The existence of CBO and their effect on the positron 
time spectra can be seen from Fig. ^, showing the Fourier 
amplitudes of residuals from fits to the data based only 
on muon decay and spin precession (Eq. ||) . The Fourier 
amplitudes at the interference frequencies w c b ,hiwa dif- 
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FIG. 2: Coherent betatron oscillations (CBO) in the g — 2 
time spectra. The Fourier spectrum was obtained from resid- 
uals from a fit based on muon decay and spin precession 
(Eq. |^) alone. The horizontal modulation was at cj c bo,h/(27r) 
= 466 kHz in the year 2000, so that the interference frequency 
i^cbo,h — w„ is numerically close to u a , as indicated. The fre- 
quency uj a is determined from fits that take CBO into account. 



fer by about 20% because of the aforementioned modu- 
lation of the observed phase. The characteristic CBO 
decay times were ~ 100 /xs. 

(3) The injection of the beam in narrow bunches into 
the storage ring resulted in a strong modulation of the ini- 
tial part of the positron time spectrum with the 149.2 ns 
cyclotron period. This effect was removed from the anal- 
yses by uniformly randomizing the recorded start time 
for each beam pulse over one cyclotron period. 

(4) Losses of muons during the data collection were 
minimized by controlled scraping [E| of the beam be- 
fore the data collection started. The time dependence of 
small, residual losses was determined by forming triple 
coincidences from scintillating detectors mounted to the 
face of three adjacent calorimeters. Corrections were 
made for accidental coincidences and proton background. 
The time spectrum of lost muons was measured, up to a 
normalization, at five locations in the storage ring. 

(5) Average calorimeter timing shifts were determined 
with a pulsed UV laser system to be typically smaller 
than 4ps over the first 200 /is of data taking. Detector 
gain changes are determined from the average observed 
positron energy as a function of time after beam injec- 
tion. ^From 50 /is on, 5 /its after the last of the calorime- 
ters was gated on in 2000, the gains of all but 2 of the 
calorimeters were stable to within 0.2% over the 10 di- 
lated muon lifetimes of data collection. 

The event sample collected in the year 2000 amounts, 
after data selection, to 4 • 10 9 positrons with energies 
greater than 2 GeV in the time region 50 /its to 600 (is 
following muon injection into the storage ring. A sta- 



tistical uncertainty in the fitted frequency u> a of about 
0.7 ppm results. To assess systematic uncertainties to 
an adequate precision, distortions of the decay positron 
spectrum in Eq. ^ have been studied in the full range 
of observed energies. In addition, the time spectrum for 
positron energies exceeding 2 GeV has been studied, as 
well as two representations of the time spectra that sup- 
press periodic distortions with frequency w c bo,h and slow 
distortions, respectively. In total, four independent and 
complete analyses of u a were performed. 

In the first analysis the data were fitted in 0.2 GeV en- 
ergy intervals in the range 1.4-3.2 GeV for each detector 
separately. This results in 198 independent fits. The 
pilcup contribution to each spectrum was determined 
from the WFD recordings using a variant of the tech- 
nique used before. For each energy interval and each 
calorimeter station, pileup was statistically constructed 
and fitted to an appropriate functional form. The re- 
sulting parameters were then incorporated into the fit- 
ting function used to describe the positron time spectra. 
Horizontal CBO were incorporated in the function fitted 
to the data through modulation of the observed number 
of counts, asymmetry, and phase. The function fitted to 
the calorimeter data in each energy interval is given by 

N(t) = N (t; E)e- t/M [1 + A(t; E) sm(uj a t + (f> a (t; E))] , 

(3) 

in which the time dependences of No, A, and <p a are all 
modulated as ~ 1 + Ai{t) sin(w c bo,h£ + 0j)j where i stands 
for any of the three separate modulations of No, A, and 
4> a - The modulation amplitudes were determined empir- 
ically from the data and follow primarily an exponential 
with a characteristic decay time of approximately 100 /is. 
Muon losses were constrained to the shape of the mea- 
sured losses. The spectrum described by Eq. ^ is multi- 
plied by a time-dependent factor which accounts for these 
losses. Fit start times for each calorimeter station and en- 
ergy interval were chosen at times after which the quality 
of fit has become constant. Slowly vanishing distortions 
of the exponential, due to imperfectly corrected energy 
scale variations, pileup, or muon losses, are thereby re- 
moved. 

In the second analysis the time spectrum of all 
positrons with energies larger than 2 GeV and arrival 
times in the region 49 fis to 600 [is after muon injection 
was considered. The pileup contribution to the spectrum 
was determined from the data and corrected for as in 
Ref. ||. The effects of horizontal CBO were incorpo- 
rated as modulations of the observed number of counts 
and asymmetry, whereas the modulation of the observed 
phase was neglected in the fitting function. The CBO 
modulation envelope was determined by partial Fourier 
integration and confirmed from the modulation of the 
observed average energies for equal phases of the (g — 2) 
oscillation. The time dilated muon lifetime 7T was fixed 
to the expected value. The treatment of muon losses fol- 
lowed the one described above. Slowly vanishing distor- 
tions of the positron spectra were parametrized empiri- 
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cally, and were included in the fitting function to improve 
the overall quality of fit. Their correlation with the fitted 
frequency Lo a is weak. 

In the third analysis, as in the second analysis, uj a was 
determined from a multiparameter fit to the summed 
time spectrum of positrons detected with all calorime- 
ters. Alternative methods were used to assess systematic 
uncertainties. In particular, to study the systematics as- 
sociated with CBO, the data were strobed at fixed phases 
of the horizontal CBO modulation. Since w c bo,h > 2xw , 
the (g — 2) signal can be fully determined from the result- 
ing spectra. Eq. | was fitted to the spectra after strobing, 
giving a result for u> a which is less sensitive by a factor of 
about three to CBO modulations of the asymmetry and 
phase than results from fits to the unstrobed spectrum. 

In the fourth analysis the ratio introduced in Eq. 6 of 
Ref. @, r(t) ~ A(E) sm(u) a t+(f> a (E)), was fitted. The ra- 
tio is largely insensitive to changes of observed counts on 
time scales larger than r a = 2ir /us a ~ 4 /us. The effects of 
asymmetry and phase modulation were not explicitly in- 
cluded in the fitting function. This analysis used summed 
calorimeter spectra, in which the asymmetry and phase 
modulation effects are suppressed by an order of mag- 
nitude as compared to summed results from fits to in- 
dividual calorimeter spectra. The suppression originates 
in the symmetry of the experiment design, and has been 
studied extensively by including the effects in the fits, by 
determining the energy modulation, and by Monte Carlo 
simulations. 

The internal consistency of the results was verified in 
various ways. For all analyses the fitted frequency uj a was 
found constant with the fit start-time, up to statistical 
fluctuations. The Fourier transform of residuals from fits 
incorporating asymmetry and phase modulation exhibits 
none of the CBO related interference structure in Fig |^. 
Fig. |a shows the results for fits to the spectra from indi- 
vidual detectors from the analysis described first. Fig. |b 
shows the results on u a versus positron energy. 

The results from the analyses are found to agree, on u a 
to within 0.4 ppm. This is within the statistical variation 
of 0.5 ppm expected from the use of slightly different data 
selections and the different treatment of the data in the 
respective analyses. Each of the analyses gives results 
with comparable total uncertainties in u> a . The analy- 
sis which incorporates asymmetry and phase modulation 
in the fitted function results in the smallest systematic 
uncertainty, primarily since CBO related systematic un- 
certainties are at the level of 0.05 ppm in this approach. 
The trade-off is a slightly enlarged statistical uncertainty, 
since more CBO related parameters are fitted. The sys- 
tematic uncertainty due to CBO is larger than our 1999 
estimate || which did not include the effects from asym- 
metry and phase modulation. Neither the 1999 value of 
dju, nor its total uncertainty were changed by this im- 
proved estimate. Other uncertainties are comparable or 
smaller. 

We combine the present results to w a /(27r) = 
229 074.11(14)(7)Hz (0.7 ppm), which includes a correc- 
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FIG. 3: a) The frequency w a /(27r) determined from fits to the 
individual calorimeter time spectra. Data from calorimeters 2 
and 20 were discarded, as in the analysis of our 1999 data, b) 
The fitted frequency uj a /{2n) versus positron energy. These 
results come from the analysis described first in the text. 



TABLE II: Systematic uncertainties for the ui a analysis. 



Source of errors 


Size [ppm] 


Coherent betatron oscillations 


0.21 


Pileup 


0.13 


Gain changes 


0.13 


Lost muons 


0.10 


Binning and fitting procedure 


0.06 


Others 1 


0.06 


Total systematic error on tu a 


0.31 



* AGS background, timing shifts, E field and vertical os- 
cillations, beam debunching/randomization. 

tion of +0.76(3) ppm for contributions to Eq. [[[caused by 
vertical oscillations and, for muons with 7 =^ 29.3, by hor- 
izontal electric fields. The stated uncertainties account 
for strong correlations among the individual results, both 
statistical and systematic. Table || lists the systematic 
uncertainties in the combined result with these correla- 
tions taken into account. 

After the uo p and uj a analyses were finalized separately 
and independently, a M was evaluated. The result is 



R 



= 11659 204(7)(5) x 10~ 10 (0.7ppm) 



(4) 



in which R = uj a /uj p and A = M M /^p = 
3.183 345 39(10) Q. This new result is in good agreement 
with previous measurements [@. || |, and reduces the 
combined uncertainty by about half. The present world 
average experimental value is 

a M (exp) = 11659 203(8) x 10~ 10 (0.7ppm), (5) 

which is driven by our determinations of a M in Refs. [| |] 
and Eq. |, and accounts for correlations between system- 
atic uncertainties. 
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Additional data on e + e~ collisions [[l2| 16 and on r- 
decay [jlTj have been published, and are being considered 
in future evaluations of a M (had, 1). New data can be 



expected from the Frascati 
factories. 



factory and from the B 



The three most recent measurements of along with 
the above standard model prediction are shown in Fig. ^|. 
The present experimental uncertainty is about half the 
size of the weak contribution to a^(SM). The difference 
of a^(exp) and a M (SM) is about 2.6 times the combined 
experimental and stated theoretical uncertainty. 



FIG. 4: Recent measurements of a M , together with the stan- 
dard model prediction using the evaluation in Ref. |J of 
a M (had, 1) from e + e~ and r decay data. 



The theoretical value of in the standard model (SM) 
is determined from a M (SM) = a M (QED) + a M (had) + 
a M (weak), in which a M (QED) = 11658 470.57(0.29) x 
10" 10 (0.025 ppm) § and o /1 (weak) = 15.1(0.4) x lO" 10 
(0.03 ppm) |0. The hadronic contribution a M (had) re- 
ceives its leading contribution from a At (had, 1), which is 
currently evaluated to be a M (had, 1) = 692(6) x 10~ 10 
(0.6 ppm) tel. The two more recent, published evalua- 
tions |h], |ll[ used recent, preliminary data from Novosi- 
birsk which have since been superseded |1| , so we choose 
not to use their values for comparison. Higher order 
contributions include a M (had,2) = -10.0(0.6) x 10" 10 
|Q and the contribution from hadronic light-by-light 
scattering, which we now take to be a^(had,lbl) = 
8.6(3.2) x 10~ 10 juj. The chiral perturbation calculation 
of Ref. [ fl5| provides an estimate with larger uncertainty. 
Hence, the value of a M (SM) is currently evaluated to be, 

a M (SM) = 11659177(7) x 10~ 10 (0.6 ppm) (6) 



In 2001 data on fj,~ were obtained. Approximately 
3 x 10 9 decay electrons were observed. Field focusing 
indices n = 0.122 and n = 0.142 were used. Measurement 
of an- will provide a sensitive test of CPT violation and 
also an improved value of a^. We plan further data- 
taking with fj,~ to obtain an additional 6 x 10 9 counts. 
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